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@ Signaltrennungsverfahren und -einrichtung fur nichtltneare Mischungen unbekannter Signale 

(57) Die Erfindung ist ein Verfahren und eine Einrichtung zur 
Separierung eines Signalgemisches, bestehend aus der 
nichtlinearen Oberiagerung von M unbekannten Signalen. 
dadurch gekennzeichnet, das N (N ^ M) unterschiedliche 
Mischungen der M Quellsignale inklusive eines eventuellen 
Storsignals einer Signalauswerteeinrichtung zugefuhrt wer- 
den, die durch eine statistische Analyse der Signale die 
nichtlinearen Ubertragungsfaktoren bestimmt und mtt die* 
sen errechneten Faktoren die Mischung ruckgangig macht, 
so da& die N Ausgange der Signaltrennungseinrichtung die 
M Quellsignale ohne Oberlagerungen enthalten. Dadurch 
wird erst ma Is die Behandlung nichtlinearer Gemische mog- 
lich, wobei nichtlinear bedeutet, daft die Quellsignale durch 
ein unbekanntes nichtlineares System gemischt werden. Das 
unbekannte System wird durch eine Volterra-Reihe be- 
■ schrieben. und die Signaltrennungseinrichtung bestimmt die 
I Koeffizienten der Volterra-Reihe. Mit derer Kenntnis ist eine 
Entmischung mdglich, auSerdem konnen die Koeffizienten in 
! weiteren Analysen zur Orts- oder Geschwindigkeitsbestim- 
mung der Signalquellen benutzt werden. 
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Beschreibung 
1. Technisches Gebiet 

5 

In der MeQtechnik besteht fast immer das Problem, daB bei einer Messung nicht nur Signale einer einzelnen 
Quelle, sondern oft nur Mischungen aus den Signalen mehrerer Quellen durch die Detektoren aufgenommen 
werden konnen. 

10 2. Stand der Technik 

Eine Signaltrennung mit mdglichst wenig a-priori- Information iiber die Signaie wurde von Herault und Jutten 
in den folgenden Veroffentlichungen behandelt: C Jutten and J. Herault, "Blind separation of sources, part I: An 
adaptive algorithm based on neuromatic architecture*, Signal Processing, vol 24, No. 1 pp 1 — 10, July 1991; 

15 P. Common, C Jutten, J. Herault, "Blind Separation of sources, part II: Problems statement". Signal Processing, 
vol 24, No. 1, pp 1 1 -20, July 1991, EL Sorouchyari, C. Jutten, J. Herault, "Blind Separation of Sources, part III: 
Stability analysis, ", Signal Processing, vol 24, No. 1, pp 21—29, July 1991. Dabei wird das Problem von mit 
konstanten Faktoren linear uberiagerten, zeitlich unverzogerten, unabhangigen Signalen beschrieben. Wie sich 
herausstellt, funktioniert das Verfahren nicht mit GauQverteilten Signalen, das sind zum Beispiel Sprache oder 

20 Musik. 

In der Veroffentlichung von Molgedey und Schuster (Molgedey L. und Schuster H. G. 1994: Separation of a 
Mixture of Independent Signals using Time Delayed Corelations, Phys. Rev. Lett. 72, 3634—3638.) wird die 
Trennung fur unkorrelierte Signale, die zeitlich unverzogert linear uberlagert sind, behandelt, die auch fur 
GauB'sche Quellen benutzt werden kann. AuBerdem wird ein Verfahren zur Trennung von Mischung skizziert, 

25 bei denen die Quellsignale (zeitlich unverzogert) multipliziert sind. 

Gilles Burell (Burell, Gilles 1992: Blind Separation of Sources: A Nonlinear Neural Algorithm, Neural Net- 
works, Vol. 5, pp 937—.) schlagt einen weiteren Ansatz fur die Trennung von Gemischen aus unabhangigen 
Signalen vor, in dem ein neuronales Netz speziell fur eine bekannte Art der nichtlinearen, zeitlich unverzogerten 
Mischungen konstruiert wird. Das Verfahren zum Finden der richtigen Parameter fur die Entmischung ist dem 

30 Backpropagationalgorithmus ahnlich. 

In dem European Patent application number 93630032.6, date of filing 10.04. 1992 by the inventors Ehud 
Weinstein and Daniel Yellin wird eine Signaltrennung fur den Fall einer durch lineare Filter entstandenen 
Mischung unabhangiger und unkorrelierter Quellsignale beschrieben, und es werden zwei Berechnungsmdglich- 
keiten fur die gesuchten Filter gegeben: 1. Eine iterative Losung zweter "standard least square" Algorithmen und 

35 2. em rekursiver Algorithraus, der im Prinzip dem ersten Equivalent ist, aber weniger Rechenkapazitat bendtigt, 
dafur aber auch nicht exakt ist Damit wird der von Piatt und Faggin in (J. C Piatt and F. Faggin: Networks for 
the Separation of Sources that are Superimposed and Delayed, ed. John E. Moody "Advances in Neural Informa- 
tion Processing" Vol.4, Morgan- Kaufman, San Mateo (1992).) vorgeschiagenen Algorithmus begriindet und 
erweitert. In keiner uns bekannten Veroffentlichung wird das Problem einer a priori unbekannten nichtlinearen 

40 Mischung (Mischung durch ein beliebiges nichtlineares System) angeschnitten oder gar gelost 

3. Zugrundeliegendes Problem 

Die Erfindung betrifft ein Verfahren sowie eine Einrichtung zur Signaltrennung eines Signalgemisches der 
45 Form 



i=l J ° 

J^. f*° f**° 

+ zZ I I ^2.Jk./( t/ i,u 2 )xjt(i - ui)xi(t - u 2 )dn x du 2 



+ 

55 M 



r m,(/ - «l)...a*m n (/ - U n )dlH ...du ny 



das aus M unbekannten Quellsignalen x ( (t) . . x M (t) besteht, die mit unbekannten Gewichtsfaktoren g 5 ^ j,,.., j n in 
der angegebenen Weise gemischt sind Ein Schema dazu steilt Fig, Nr. 1 dar. 1 steht fur den unbekannten^Teil, 
bestehend aus den Quellprozessen, deren Signale durch ein nichtlineares System G 2 gemischt werden. 3 zeigt 
den bekannten Bereich an. In 4 werden aus einer Analyse der MeBdaten die Quellsignale entmischt und stehen 
65 als Ausgange u, und u 2 zur Verfugung. Das Verfahren kommt fast ohne Annahme iiber die Quellsignale und das 
nichtlineare Mischungssystem aus, bis auf die Bedingung, daB eine statistische Analyse der empfangenen Mi- 
schungen aus den Quellsignalen y^t) geniigend Information fur eine Trennung Hefert. Beispielsweise wird die 
zeitlich unverzogerfe Oberlagerung zweier GauBverteilter weiBer Rauschsignale ni(t) und n^t) mit unterschied- 
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'I 

lichen konstanten Faktoren C12 und c 2 » in der Form 

yi(t) = ni(t)+ci 2 n2(t) 
y2 (t) = n2<t) + c ( 2ni(t) (1) 

5 

nicht moglich sein, da (yiWy^t + x)) = 0 fur t^O und alle Korrelationen der Form (yVOy^t + x)) = 0 fur 
m + n>2 

Die Beschrankungen werden im Verlauf dieser Beschreibung klarer werden, aber da eine Veranderung der 
g s n,mi,..,m n durch eine geanderte Aufnahmetechnik moglich ist und da man es in den meisten Fallen mit der 
Trennung von Signalen, die Information enthalten, d. h. bei denen eine statistische Analyse hoherer Ordnung io 
Ergebnisse liefert, wie z. B. eine Trennung von Mischungen aus Sprache und Musik, zu tun hat, werden die 
Beschrankungen meistens ausraumbar bzw. vernachlassigbar sein. 

4. Erfindnung 

15 

Eine beliebige nichdineare Mischung kann durch die folgende Form beschrieben werden: 
M _ 

20 



//.(') = £ jf" <jij(i'i)rjlf -u x )<lu x 

roo r<x> 



25 
-r . . - 

+ H L •••/ ^l, 

m 1 ...m l( =:] ^0 JO 

x mt x mn (f - u n )rf«i . - . du n , (2) * 30 

wobei die M Queilsignale Xt(t) . . . xm(0 unbekannt und mit unbekannten Gewichtsfaktoren g 5 n,ji t ..,j 0 in der 
angegebenen Weise gemischt sind und wobei teR ist Dies ist die Form einer Volterra-Reihe. Entsprechend 
einem Theorem von Frechet kann jedes stetige Funktional rait einer beliebigen Genauigkeit durch eine Volter- 35 
ra-Reihe approximiert werden. Zu dieser Reihe gibt es auch ein zeitdiskontinuierliches Analogon, das eine 
Anwendung der Erfindung auf digitate Systeme zulaBt. 

Im folgenden wird nur der Fall betrachtet, in dem xj(t), y$(t), t € R gilt Eine Formulierung mit hoherdimensio- 
nalen Variablen ist moglich, aber der eindimensionaJe Fall enthalt alle wesentlichen Verfahrensschritte. 

Das in dieser Erfindung beschriebene Verfahren besteht im wesentlichen aus zwei Schritten: Erstens werden 40 
die durch den wahlbaren Grad der Nichtlinearitat bei der Mischung eindeutig bestimmten nichtlinearen Glei- 
chungen fur ein gleitendes Zeitfenster gelost und die Losungen werden iiber die Zeit gemittelt, und zweitens 
wird ein aus einer genugend groBen Anzahl von unterschiedlichen Kumulanten der geschatzten Ausgangssigna- 
le ui(t) gebildetes Potential minimiert, wobei die zur Berechnung des Potentials notigen Werte aus einem 
gleitenden Zeitfenster einer wahlbaren Lange stammen. 45 

Zu erstens: 

Eine mehrdimensionale Fouriertransformation von 2 ist moglich, und man erhalt: 



W) = I! )-Y j (u; l ) + .- ■ 
i=i 

M 

+ Y. G «..n, u, n )A' m ,(u;,)...A' m „(u.„). 

1711 »< tl = l 

Genauso ist eine zeitvariable Fouriertransformation mdglich, in der jeweils ein zeitlich begrenzter Ausschnitt 
der GroBe TeR aus dem Gesamtverlauf der Signale herausgeschnitten und fouriertransformiert wird. Dann 
erhalt man: 

VA^T) = £<7J J (wi)A' i (^.,7 , ) + ... 

A/ 

(W, , . . . , UfcJA'm, (W, , T) . . . X m „ (W n , T). 
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Dabei ist vorausgesetzt, daB sich das System G so langsam andert, daB diese Anderung bei der Berechnung der 
gesuchten Mischungsfaktorcn vernachlassigt werden kann. 

Urn die Bedingungen zu tormulieren, die die Quellsignale erfiillen mussen, damit das Verfahren im ersten 
Schritt angewendet werden kann, werden nun folgende Variablen definiert: 

5 ijJcNftNuNaeN 

Ti jk eR 

(ii,...iN 0 )e{peN|l<p<M| N o 
I0 (j.,...jN 1 )eN N » 
(ki,..,k N2 )eN N 2 

(©i,jl,...aHN<tfN|>6R N 0 N l 

( TijJtk,, . . , TiN0jNlkN 2 ) € R N 0NlN 2 

is Die Quellensignale mussen die foigenden Bedingungen erfiillen fur einen genugend groBen Bereich der oben 
definierten Variablen: 



ctiw( (3) 

A», {uixjfif! » T — n, y #Vl a, ),-... A„ (w.-jj^ , T — ^ ) t 
1 = 0. 

35 auBer falls 

i t = i 2 . . . = j No 
und 

6>ilJl = . . . = OiNOfjNl 
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gleichzeitig gilt 

Nun geht es darum, die in der obigen Gleichung stehenden unbekannten GroBen zu berechnen. Man benotigt 
45 dieselbe Anzahl an Gleichungen wie an unbekannte GroBen. Daher soil mit obiger Gleichung fur ausreichend 
unterschiedlich viele Verzogerungszeiten x 0 bis te(E€N) ein Gleichungssystem aufgestellt werden, d h.; 



(KoC^rjy^^r + To)) = 



(Y S0 (^T)Y Si (Z,T + t e )) = ... (4) 

55 oder entsprechende Gleichungen, die sich durch das Bilden der Erwartungswerte von beliebigen Produkten der 
zeitvariablen Fouriertransformierten der empfangenen Signale Ys(g>i, . . , o n , T) und der zeitverzogerten zeitva- 
riablen Fouriertransformierten der empfangenen Signale Y s (o>t, . . o>n, T— Xj) ergeben. Dadurch erhalt man ein 
nichtlineares Gleichungssystem, das numerisch gelost werden kann. Naturlich kann man durch die Beriicksichti- 
gung weiterer Werte fur die Verzogerungszeiten xj mehr Gleichungen als unbedingt notig aufstellen und dann 
6o die uberbestimmten Losungen einer Mittelwertbildung unterziehen. 

An einem Beispiel soil nun die Verfahrensweise dieses ersten Schrittes gezeigt werden. 
Das betrachtete System sei derart nichtlinear, daB gilt: 



65 
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G!. a («,)X,(u; l .r) 

ci. 2 . 1 (u,. i .- 2 )X.( t ' 1 .r)^.(^- F): 
y J («_- l ,ej 3 .7") = x 2 {<-i.T) -5- 

6^,(^)^,(0., : n 

C\ aA («;, , to, )X 2 («i , r)Xi 7*). 

Dann ergeben sich durch das Bilden von Kremkorrelationen die folgenden Gleichungen (0<j<E) fur 
©2€Rund©i9fc©2: 

<Yi(©,,© 2 ,T)Y 2 (©i,© 2 ,T+xj)) = 
G»i.i(«fliXXi(B>i,1)Xi(a>i,T+T^ + 
G' i j(a>i)(X2(a>i,T)X3(<tti, T+ ij)) + 

G'zUo),, co 2 )G 2 2.i5(o) l , a>2)(Xi(o>i, T)X,(©,,T + Tj)XX2(«B 2 ,T)X2(ojj, T+ xj)) + 
G'2Ai(o>,, © 2 )G 2 2Ai(coi, o> 2 XX 2 (a>,, T)X2(©,,T+xj)XX,(© 2 ,T)X,(© 2 , T+ xj)); 

(Y,(©,,© 2 ,T)Yi(©,,©2,T+Xj)) = 

(X,(©,,T)X I (©,,T + x)) + 

(G« t^(<o))^X2(a>i, 1>X2(«>i. T+ xj)) + 

(G'2.1^{<0, f ©z)) 2 ^©,, T)X,(©,, T + TjJXX^CO* T + Xj)> + 

(G 1 2Ai(©,,to 2 )) 2 ^2(o>i.T)X 2 (<o 1 ,T+T J )XX,(ffl 2 ,T)X,(a> 2 ,T+x j )>; 

(Y^co,, 0)2, "HY^©, , ©2, T + xj)) = 
(X^©,, T)X2(<i)i , T + Tj)) + 
(tfu^TOK T)X,(a>,, T+ xj)) + 

(G 2 2,U&)i, ©z)) 2 ^©,, T)X,(©,, T + TjJXX^©* T)X2(© 2 , T + xj)) + 
(G^wKffli. ©z)) 2 ^©,, TJXzCo),, T + Tj)XX,(©2. TW©* T + xj)). 

Entsprechend ergeben sich 3 weitere Gleichungen, wenn ©, und © 2 vertauscht werden. 

(Y,(© 2 .©i,T)Y 2 (©i,© 2 .T+Tj)) = 
G 2 i,i(© 2 )(X 1 (© 2 , T)X,(© 2 , T + xj))+ 
G" i^(<D2XX2{<a2. "DX^a* T + xj)) + 

G'z,^©* tsuy&v d*** o>iXXi(a)2, T)X,(© 2 , T + XjJXX^©,, TJX^©,, T + xj)) + 
G»2Ai(<»2. ©iJG^iC©* ©,XX2(© 2 , T)X2(©2. T + Xj)XX,(©,,T)X,(©,, T+ Xj)); 

(Y ,(© 2 , ©,, T)y,(© 2 , ©,, T+ xj)) = 

(X,(© 2 . T)X,(© 2 , T + xj)) + 

(G 1 ijfo^ftXrfcta, TJXK©* T + Tj)) + 

(G' 2 ,iX© 2 .©i)) 2 (X l (©2,T)X,(©2,T+xj)XX 2 (© l ,T)X2(©,,T+Tj)) + 
(G'2Ai(o> 2 . ©i)) 2 ^^, T)X 2 (©2,T+Xj)XX 1 (©,,T)X,(©, > T+ Xj)); 

(Y 2 (©2,©i,T)Y 2 (© 2 ,©,,T+xj)) = 
(X2(© 2 , TJXzf©* T + Xj)) + 
(G 2 ,,i(© 2 )) 2 (X,(© 2 , T)X,(© 2 , T + Xj)) + 

(G 2 2,U© 2 , ©,)) 2 <X I (© 2 , T)X,(© 2 , T + xj)XX2<©, , T)X2{©,, T + xj)) + 
(G 2 2Ai(© 2 , ©i)) 2 ^©* T)X2(© 2 . T +Xj)XX,(©,, T)X,(©,, T+ xj)). 

Fur ein Paar (©,, © 2 ) mit ©i <t © 2 ergeben sich 6E Gleichungen fur 12 + 4E unbekannte GroBen. 
Im Fall, daB ©,, © 2 SR mit ©, =© 2 ist, ergeben sich durch das Bilden von Kreuzkorrelationen die folgenden 
Gleichungen (0 < j < E): 

(Y,(©,,©,.T)Y2(©i,©i,T + Xj)> = 
(G 2 U (© 1 XX,(©,,T)X,(© I ,T+Xj)) + 
(G' i 2 (© l XX 2 (© 1 ,T)X2(©,,T+ xj))+ 

(G'2, I2 (©i» ©iJG 2 *,;*©,, »i)[Xi(tti, T)X,(©,. T+XjJ^X^©,, T)X^i. T+ Tj))+ 
G l i2 ,i(^i.<a,)G 2 2Ai(©.,©iXX2(©i.T)X2(©,,T + Xj)XX 1 (©,,T)X,(©i,T+xj)) + 
G I i i J (©,,©,)G 2 2ii(co!,© 2 XX 2 (©i,T)X2(©,,T+xj)XX,(©,,T)X,(©,,T+Xj))+ 
G l i2,i(©i. ©,)G 2 zi j(©i, ©,XXj(©i, T)X 2 (©,, T +xj)XX,(©,, T)X,(©,. T+ Xj)); 
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(Yifoi.OuTJYiCoi.cowT+Tj)) = 
(X|(<oi,T)X((a)| f T + T j )> + 
(G 1 Uo.TOa),, TJXKo, , T + xj)> + 

5 (G i 2A i(to, f coi))pC 2 (couT)X 2 (a)i ( T+T j )XX I (a)i,T)X 1 (co,.T+T i )> + 
G'^ai.oiJG'uito.ah^ 

(Y 2 (<Di,(i>i,T)Y2(cDi f 0>i,T+Tj)) = 

io (X^co i , T)X 2 (a)u T + xj)) + 

(G 2 ij(ct>i)XXi(a>i,T)Xi(coi f T + Tj)) + 

(G^tti, coi))2(X t (a> lr TPC(<o lf T+ xjJXX^co,, T)X 2 (a>,, T+ Xj)) + 
(G^Koi. ©i)fpC2(a), f TJXKa),, T+ Tj)XX,(a>i f T)X,(to ( , T+xj)) + 
G^^ai.aOG^ttui,^ 
15 G^^ai.aO&ui^i^ 

Fur ein Paar (g>i, tt> 2 ) mit o>i = o 2 ergeben sich 3E Gleichungen fur 6 + 2E unbekannte GroBen. 
Beim Aufstellen dieses Gleichungssystems wurde davon Gebrauch gemacht, daB die Quellsignale die folgen- 
den Bedingungen entsprechend Gleichung (3) erf ullen: 
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cum(Xi (wi , T), X 2 (u> 7 , T — 






0 


cttm(A' I (a Jl ,r),.Y|(w„r- 


t),X 2 {uj 2 ,T - 


*-)) 




0 


am(A'j(«tfi , 7*), A" t (w, , T - 


t),X 2 (uj 2 ,T — 


*-)) 




0 


cum( Xi («, , T) , X 2 (ui , T - 


r),-Yi(wi,r- 


r)) 




0 


cum(A' a (w, , T), A 2 (u>, , T - 


r),X x (u.- 2 .T- 


r)) 




0 


cum(A' t (u,, , T), A'ato, T), .V,(u, , r - 


t),X 2 (uj 2 ,T - 


r)) 




0 


cum( A,(w, , T), A' 2 («a, T), X 2 (u u T - 




r)) 




0 


, T), A', (a*, T), A' 2 (wi , T - 


r),A'i("2,r- 


r)). 




0 


cum(A' 2 (u» 1> D, A"i(a.«2, T), X^ U T - 


T),A' 2 («^,r- 


')) 




0 



fur alle o>i, o> 2f te R, und 

40 cum (X,(coi, T), X,(a> 2 , T-x)) = 0 
cum (X^o,, T), XKcoj, T-x)) = 0 

fur alle Oi o> 2l x € R mit cdi ^= o> 2 . 
Damitgiltdann: 

45 

(K t ((a u T)X2(o) 2t T)Xi(o>i , T— x)X^ 2t T-x)) = 
(X,(6) If T)X ( (o> If T- x)XX2(<o 2f T)X2(a> 2> T- x)) ; 
(X^oi, T)Xi(o 2f T)X2(oi, T-x)X,(a) 2 , T-x)) = 
(X^o,, T)X2<<d 1( T-x)XXi(a> 2 , T)X,(a> 2 , T- x)>; 
so (XiCtou-nX^coi^X^o^T-xJX^OuT-x)) - 
(X 2 ((o l ,T)X 2 (a)»,T-x)XX l (co,,T)X I (a>uT-x)); 

fiir alle o>i, o> 2 € R- AuBerdem gilt: 

55 {Xi(o>i f T)X 2 (o> 2> T)X 2 (o3|,T-x)X,(6) lf T-x)) = 0; 
(X^o^TJXiCco^T^XKcoi.T-xJX^cowT-x)) = 0; 

fiir alle o>t, o> 2 e R mit cot ^ a> 2 . 

Aufgrund von fehlerhaften Aufnahmen (zusatzliches Rauschen) bzw. Fensterungseffekten durch die zeitvari- 

eo ablen Fouriertransformation muB diese Berechnung fur ein in geeignet groBen Schritten uber die Sign ale 
bewegtes Fenster wiederholt werden. Die entsprechenden Ergebnisse werden dann gemittelt, und die gemittel- 
ten Werte nahern sich so den wahren Werten an. Diese Annaherung kann aber viele Daten bendtigen. Eine 
groBe Menge an Daten kann zu Problemen bei einer Realisierung fuhren, daher wird diesem ersten Schritt ein 
zweiter nachgesteilt 

65 Zu zweitens: 

Dieser zweite Schritt ist eine Anlehnung an die Methoden, die in der Signalverarbeitung mit neuronalen Netzen 
angewendet werden. Man konstruiert sich eine Kostenfunktion, die minimiert wird. Wenn das globale Minimum 
erreicht ist, hat man die optimalen Werte, in diesem Fall der Obertragungsfaktoren, gefunden. Fur die Minimie- 
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rung einer Kostenfunktion stehen beim heutigen Stand der Technik eine groBe Anzahi von Methoden zur 
Verfugung, z. B. Gradientenabstieg oder Monte-Cario-Simulationen. Ein Nachteil aller Verfahren ist aber, daB 
sie in lokalen Minima der Kostenfunktion steckenbleiben kdnnen; der Vorteil unserer Erfindung ist, daB durch 
die Startwerte aus dem ersten Schritt die Minimierung schon in der Nahe des giobalen Minimums beginnen 
kann, denn bei dem Ende der Berechnungen des ersten Schrittes liegen fur die Obertragungsfaktoren geschatzte 5 
Werte vor, die sich nur wenig von den gesuchten Werten fur die wahren Obertragungsfaktoren unterscheiden 
und die fur eine bessere Trennung weiter optimiert werden. 

Mit der Methode der Potentialminimierung werden Effekte von zusatzlichen Rauchprozessen verringert 
Urn die Bedingungen zu formulieren, die die Quellsignale erfiillen mtissen, damit das Verfahren im zweiten 
Schritt angewendet werden kann, werden nun folgende Variablen definiert: io 

iJ,No,NieN (5) 

(i,...JN 0 )e{peN|l<p<M} N o 

Oi,...jN,)€N N l I5 

(T il i 1 ....T i NojNi)eR N ° Nl 

Die Quellensignale raiissen die folgenden Bedingungen erfullen fur einen genugend groBen Bereich der oben 
definierten Variablen: 

20 

cum( (6) 

x ii ~~ r t'i Ji )i - - - ) * r ti ~ T ii Jjvj ), 
• r t,v 0 ~~ T iN 0 Ji )»■••-» X irr Q ~~ T tv 0 Jjv, ) 

) = o, 

auBer falls 
ii = 12 = ... = iN 0 

gilt. 35 

Fur den hier betrachteten Fall, die richtigen Obertragungsfaktoren zu finden, bilden wir eine Kostenfunktion 
derart, daB die Ausgange eines adaptiven Filters u;(t)dieselben Bedingungen erfullen rniissen wie die Quellsigna- 
le xi(t), & h. far eine genOgend groBe AnzahJ von Werten der in 5 definierten Variablen werden die Quadrate der 
Kumuianten aus 6, dabei werden die Quellsignale x durch die Ausgange des Filters u ersetzt, aufsummiert, so daB 
eine Kostenfunktion entsteht, die minimiert wird. Optimiert werden dabei die Werte der Obertragungsfaktoren. 40 
Nach jedem Optimierungsschritt rniissen die Werte fiir die Ui(t) neu berechnet werden. Diese Berechnung kann 
im Zeit- oder Frequenzbereich erfolgen. 

Zum Beispiel werden fiir den Fall, in dem folgende Mischung vorliegt: 

Y,(o>) - X^caJ + d^X^cojH- K(co)X,(co)X2(o); 45 
Y 2 (o>) = X2(ti>) + G 2 i(o>)X|(o))+K(co)X,(G>)X2(a)) 

die Ausgange mit den geschatzten Werten fiir die Obertragungsfaktoren G*i2, G*2i und K* berechnet durch 
Auflosen des Gleichungssysteras: 

Y,(o) = U^oJ+G'^coJU^ + K^coJUU^UKti)); 50 
Y2<a>) = U 2 (a)) + G*2i(o)Ut(co) + K*(ci>)U 1 (co)U2{a)) 

nach Ui(o) und U2(o>). 

Die durch das adaptive Filter zu minimierende Kostenfunktion V kann mit den Ausgangen Ui(t) und U2O) 55 
gebildet werden: 

V=Y. I«tm(t£i(0,u a (« + r))l a + 

r€B 60 

(cum(ui(t) ia ,(i + r),u 2 (t + r))] 2 + 
[citm(u 7 (t), u 2 (t + r), Wl (t + t))] 2 + 

[cuni(u l (t),u 7 (t),u l (t+T),u 2 (t + T))] 2 + 65 
[cum(u 2 (0.«i(0."a(t + r),ui(t + r))j 2 , 
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wobei B eine endliche Menge reeller Zahlen ist 

5. Gewerbliche Anwendbarkeit 

5 Die durchzufuhrenden Berechnungen bendtigen eine hohe Rechenkapazitat, die aber schon heme mit Mikro- 
prozessoren bereitsteht Eine in parallelen Schritten ausgefuhrte Verwendung des angegebene Verfahrens 
eroffnet auch fur hdherfrequente Probleme ein Echtzeitrealisierung. Im offline Betrieb laBt sich das Verfahren 
auf frei program mi erbaren Mikroprozessoren implementieren. 

Anwendungsmoglichkeiten liegen in einer Signaltrennung als Vorverarbeitungsschritt fur eine automatische 
to Spracherkennung, bei der Konstruktion technischer Horhilfen, sowie in alien anderen akustischen Messungen, 
einschlieBIich die der Unterwasserakustik. 

AuBerdem sind medizinische Anwendungsmoglichkeiten bei Mehrkanalmessungen wie EEG oder MEG 
(Fig, Nr. 2) moglich, sowie bei Antennenfeldern fur die Aufnahme elektromagnetischer StrahJung. Wie in der 
Fig. Nr. 2 dargestellt werden magnetische Dipole im Gehirn 2 durch SQUIDs 3 gemessen, und die Signale 
is werden einer Auswerteeinheit 1 zugefuhrt. 

Diesen und alien weiteren denkbaren Anwendungsmoglichkeiten ist gemeinsam, daB sie der bisher linearen 
Theorie verschlossen bleiben, weil eine lineare Signaltrennung zu unbefriedigenden Ergebnissen fiihrt, denn 
aufgrund der Detektortechniken haben fast alle Empfanger eine nichtlineare Kennlinie. Weiterhin sind bei einer 
Formulierung mit mehrdimensionalen Variablen Anwendungen in der Optik bzw. Btldverarbeitung moglich, in 
20 der auch bei der Aufnahme nichtlineare Effekte auftreten. 

Daruber hinaus konnen die berechneten Werte auch zu einer Analyse des mischenden Systems benutzt 
werden, z. B. zur Quellenlokalisierung oder Bewegungsdetektion. 

6. Vorteilhafte Auswirkungen 

25 

Neben den schon in dem Abschnitt Gewerbliche Anwendbarkeit genannten moglichen Erneuerungen gibt es 
noch eine weitere Verbesserung: Urn mit der bisher linearen Methode zu arbeiten, werden Empfindlichkeitsbe- 
reiche der Detektoren auf ihren linearen Bereich verkleinert Dies ist bei der nichtlinearen Methode nicht notig. 

30 7. Ausfuhrung 

Stellvertretend fur alle anderen Probleme wird der Fall der Trennung akustischer Signale behandelt Entspre- 
chend Fig. Nr. 1 stellt der mit 1 bezeichnete Block dann einen unbekannten ProzeB bestehend aus zwei Schall- 
quellen dar, deren Signale durch den Raum zu zwei Mikrofonen ge lan gen, deren aufgenommene Signale 

35 aufgezeichnet werden und die Grundlage der der Erfindung entsprechenden Berechnungen bilden. Der Block, 
der mit 3 gekennzeichnet ist, enthalt alle bekannten GroBen. Das System G ist nichtlinear, weil zum einen die 
Raumubertragung nichtlinear sein kann, zum anderen die Mikrofone mit ihren Verstarkern eine nichtlineare 
Kennlinie besitzen. Die berechneten Signale ui und u 2 sind den Quellsignalen k x und *t ahnlich, so daB ein 
automatisches Spracherkennungssystem mit der aus der Mischung herausgetrennten Sprache funktioniert Dies 

40 ist in der Fig. Nr. 3 sichtbar. Dort sind in den oberen beiden Kasten die Mischsignale von Sprache und Musik 
aufgezeigt wie sie sich bei einer Mikrofonaufnahme in einem kJeinen Raum ergeben- Im untersten Kasten ist das 
aus dem Gemisch herausgetrennte Sprachsignal am Ausgang der Trennungseinrichtung, dessen Berechnung auf 
einem der Erfindung entsprechenden Verfahren basiert, dargestellt; mit dem untersten Signal ist eine fehlerlose 
automatische Erkennung mit einem koramerziellen automatischen Spracherkenner moglich, wahrend derselbe 

45 Erkenner bei den oberen beiden Signalen jeweils drei Fehler liefert 

8. Figurenbeschreibung 

Fig. Nr. 1: 1 steht fur den unbekannten Bereich, bestehend aus den Quellprozessen, deren Signale durch ein 
50 nichtlineares System G 2 gemischt werden. 3 zeigt den bekannten Bereich an. In 4 werden aus einer Analyse der 
MeBdaten die Quellsignaie entmischt und stehen als Ausgange ui und U2 zur Verfugung. 

Fig. Nr. 2: Magnetische Dipole im Gehirn 2 werden durch SQUIDs 3 gemessene, und die Signale werden einer 
Auswerteeinheit 1 zugefuhrt 

Fig. Nr. 3: In den oberen beiden Kasten sind die Mischsignale von Sprache und Musik aufgezeigt, wie sie sich 
55 bei einer Mikrofonaufnahme in einem kleinen Raum ergeben- Im untersten Kasten ist das aus dem Gemisch 
herausgetrennte Sprachsignal am Ausgang der Trennungseinrichtung, dessen Berechnung auf einem der Erfin- 
dung entsprechenden Verfahren basiert, dargestellt; mit dem untersten Signal ist eine fehlerlose automatische 
Erkennung mit einem kommerziellen automatischen Spracherkenner moglich, wahrend derselbe Erkenner bei 
den oberen beiden Signalen jeweils drei Fehler liefert 

60 

Patentanspriiche 
1. Verfahren zur Separierung eines Signalgemisches der Form 
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Vs(t) 




5 
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das aus M unbekannten Quellsignalen Xi(t)... XM(t) besteht, die mit unbekannten Gewichtsfaktoren is 
gV j 1? ..,j n in der angegebenen Weise gemischt sind, fur die Variable t soli geiten: teR, dadurch gekenn- 
zeichnet, daB ersten die durch den Grad der Nichtlinearitat bei der Mischung eindeutig bestimmten 
nichtJinearen Gleichungen, die durch das Bilden von einer geniigend groBen Anzahl von unterschiedlichen 
Erwartungswerten von Produkten der zeitvariablen Fouriertransformierten der empfangenen Signale 
Y s (o>i f ..^ co n , T) und der urn Tj zeitverzogerten zeitvariablen Fouriertransformierten der empfangenen 20 
Signale Y s (wt, . . , o) n , T— Tj) entstehen, fur ein gleitendes Zeitfenster geldst werden und die Losungen liber 
die Zeit gemittelt werden und daB zweitens ein aus einer geniigend groBen Anzahl von unterschiedlichen 
Kumulanten aus den geschatzten Ausgangssignalen Ui(t) und aus den urn Tj zeitverzogerten geschatzten 
Ausgangssignalen ui(t— Tj) gebildetes Potential minimiert wird, wobei die zur Berechnung des Potentials 
notigen Werte aus einem gleitenden Zeitfenster einer wahlbaren GroBe stammen. 25 

2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daB uber die Werte der Nichtlinearitat g irgendwel- 
che a-priori-Informationen bestehen, so daB die Berechnungen vereinfacht werden konnen. 

3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB M = 2 ist 

4. Verfahren nach Anspruch 1 bis 3, dadurch gekennzeichnet, daB uber die Werte der Quellsignale x 
irgendwelche a-priori-Informationen bestehen, so daB die Berechnungen vereinfacht werden konnen. 30 

5. Einrichtungen, dadurch gekennzeichnet, daB sie die Verfahren der Anspriiche 1 bis 4 in digitaler oder 
analoger Weise fur eine Signaltrennung benutzen. 
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